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1.  Introduction 

Modern  large  buildings  in  the  Tokyo  metropoli- 
tan area  usually  have  a centralized  hot-water  sup- 
ply operating  24  h a day.  Copper  pipe  is  widely 
used  in  such  systems  because  of  its  relatively  high 
resistance  to  corrosion,  coupled  with  additional 
pragmatic  merits:  it  is  easy  to  work  with,  easy  to 
install,  and  relatively  cheap. 

The  seriousness  of  Type-II  pitting  corrosion, 
however,  has  increasingly  received  high  recognition 
in  such  hot-water  supply  systems  [1,2,3].  The  need 
to  obtain  information  regarding  the  degree  of  pit- 
ting corrosion  has  increased  over  the  last  decade 
because  considerable  pipe  damage  may  require 
maintenance,  and  even  replacement,  and  in  that 
case,  proper  life  prediction  is  essential  to  pass  rea- 
sonable engineering  judgements  and  thereby  to 
perform  proper  maintenance. 

The  life  prediction  of  such  copper  plumbing  tub- 
ing can  first  be  performed  by  coupling  adequate 
nondestructive  and/or,  though  less  favorable,  de- 
structive inspection  techniques  with  reliable  statis- 
tical analysis.  It  has  been  shown  that  the  most 


promising  statistical  analysis  methods  for  such  a 
purpose  include  extreme  value  statistical  analysis 
[3,4]. 

Although  such  tools  have  become  widespread,  a 
general  method  for  evaluating  the  localized  corro- 
sion propensity  on  existing  engineering  structures 
from  limited  inspection  data,  and  concrete  crite- 
rion for  the  number  and  size  of  samples  required 
to  obtain  a reasonable  extreme  value  prediction  is 
still  not  available. 

In  the  present  study,  a set  of  pitting  corrosion 
depth  data  obtained  from  7 year  old  copper  plumb- 
ing pipe,  one  third  of  which  was  removed  from  a 
centralized  hot-water  supply  system,  was  examined 
by  extreme  value  statistics.  Emphasis  was  placed  on 
the  effect  of  the  total  number  and  size  of  maximum 
pit  depths  on  the  accuracy  of  pit  growth  prediction. 
A concept  was  ultimately  proposed  to  obtain  a rea- 
sonable prediction  of  maximum  pit  depths  while 
minimizing  the  total  sampling  area  (or,  length)  for 
analysis. 
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2.  Test  Procedures 

2.1  Test  Specimens 

Copper  pipes  totaling  8.88  m in  length  were  re- 
moved from  various  parts  of  the  centralized  hot- 
water  supply  system  in  an  11-story  multi-family 
dwelling  in  Tokyo.  The  system  consisted  of  one 
stainless  steel  storage  tank  (1.8  m^)  and  copper 
plumbing  pipe  having  an  overall  length  of  about  28 
m.  The  plumbing  material  was  JIS  C 1220T  type 
25AM  (outside  diameter  28.58  mm  and  wall  thick- 
ness 0.89  mm)  copper  pipe  for  building  use.  The 
system  had  been  operating  for  about  2600  d before 
test  piping  was  removed.  In  the  system,  water  at  a 
nominal  temperature  of  60  °C  circulated  con- 
stantly. Average  flow  (i.e.,  hot-water  consumption) 
was  around  8 mVd,  the  storage  tank  being  supplied 
automatically  with  tap  water. 

The  copper  piping  removed  was  cut  into  parts 
100  mm  long,  which  were  then  cut  in  half  to  give 
half-ring  specimens.  Each  half-ring  specimen  was 
then  completely  cleaned  ultrasonically  in  dilute  sul- 
furic acid,  followed  by  marking-off  to  divide  it  into 
10  virtual  half-ring  specimens  of  10  mm  unit  length, 
after  which  the  pit  depths  were  measured  by  using 
an  optical  microscope  of  1 |xm  precision.  By  cou- 
pling two  opposite  virtual  half-ring  specimens,  a 10 
mm  long  full  ring  specimen  was  reassembled,  and 
the  area  thus  surveyed  should  be  representative  of 
the  pipe  at  that  particular  location.  Then,  by  taking 
several  adjoining  full-rings  of  10  mm  length,  for 
each  unit  sample  sizes  of  maximum  pit  depths 
(hereinafter,  s:  unit  length)  in  the  interval  20  mm 
up  to  200  mm  were  obtained. 

2.2  Extreme  Value  Statistical  Analysis 

The  extreme  value  statistical  analysis  was  per- 
formed by  using  a commercial  available  personal 
computer  software  package,  EVANS  [5]. 

The  basic  concept  of  the  present  extreme  value 
analysis  is  briefly  reviewed  in  the  following  sec- 
tions. 

2.2.1  Extreme  Value  Probability  Plots  The 
first  step  of  the  extreme  value  analysis  included  the 
preparation  of  extreme  value  probability  plots,  that 
is,  plots  of  maximum  pit  depth  data  on  extreme 
value  probability  paper  of  the  cumulative  relative 
frequency  (F(y))  vs  maximum  pit  depth  (;c).  Maxi- 
mum pit  depths  data  were  arranged  in  order  from 
largest  to  smallest  and  assigned  a rank  number. 
The  vertical  plotting  position  F(y)  for  each  pit 
depth  value  was  calculated  by  the  averaged  rank 
method  as  follows: 


E(y)  = l-//(A  + 1)  (1) 

where  i = rank  number, 

N = total  number. 

Several  sets  of  plots  were  obtained  depending  on 
the  combinations  of  N and  s.  In  an  exemplifying 
case  of  5 = 100  mm  and  N = 1,  a total  of  10  sets  of 
plots  was  subjected  to  regression  analysis. 

2.2.2  Regression  for  the  Best  Fit  Line  Re- 
gression analysis  based  on  the  MVLUE  (Minimum 
Variance  Linear  Unbiased  Estimator)  method  was 
made  for  each  data  set  to  determine  a straight  line 
of  best  fit  to  the  plotted  extreme  values.  The  equa- 
tion of  a straight  line  is  given  as  follows: 

x — k-\-c^  (2) 

where  x = expected  maximum  pit  depth, 
y = standardized  variable, 

A = location  parameter, 
a = scale  parameter. 

The  straight  line  is  drawn  in  Fig.  1 on  extreme 
value  probability  paper.  In  the  present  analysis, 
emphasis  was  placed  on  the  optimum  combinations 
of  N and  s for  obtaining  a reasonable  estimate  of 
the  extreme  value.  Such  an  analysis  becomes  feasi- 
ble through  a thorough  investigation  of  all  the  sam- 
ples where  the  depth  of  the  actually  detected 
deepest  pit  could  be  regarded  as  the  probable  max- 
imum pit  depth,  that  is,  the  extreme  value. 


Fig.  1.  Extreme  value  probability  paper  of  double-exponential 
distribution. 
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3.  Results  and  Discussion 

3.1  Distribution  of  Pit  Depths 

Figure  2 is  a histogram  showing  the  relationship 
between  the  pit  depths  and  the  total  frequency.  In 
the  figure,  a total  of  970  depth  data  for  all  pits 
found  in  one  10  mm  long  full-ring,  removed  from 
ihe  7th  floor  in  the  building,  was  grouped  over  the 
pit  depth  ranges,  0 mm-0.019  mm,  0.020  mm-0.039 
mm,  etc.  As  is  evident  from  the  figure,  the  shape  of 
the  pit  depth  distribution  is  a bell-shaped  curve 
starting  at  zero,  rising  to  a maximum  at  around 
0.05  mm  and  thereafter  decreasing  rapidly  with  in- 
creasing pit  depths.  The  lack  of  “J”-shaped  portion 
bending  to  the  right  in  the  pit  depth  range  from  0 
mm  to  0.02  mm,  together  with  the  trailing  extreme 
portion  of  the  tail  of  the  curve  up  to  0.34  mm,  indi- 
cated that  most  of  the  small  pits  had  already  had 
ceased  to  grow  while  only  a small  number  of 
deeper  pits  continued  to  grow. 


Pit  depth  (mm) 


Fig.  2.  Pit  depths  distribution  histogram  for  10  mm  long  copper 
pipe  used  for  2600  d in  a hot-water  supply  system  operated  at 
60  °C. 

In  Fig.  3,  all  the  pit  depth  data  represented  in 
Fig.  2 were  plotted  on  a logarithmic-normal  distri- 
bution diagram.  The  apparent  linearity  of  the  plot 
indicates  that  this  distribution  applies  to  the  loga- 
rithmic-normal; and  hence  the  maximum  pit  depth 
data  obtained  in  the  present  pit  depth  survey  were 
sampled  from  the  parent  populations  with  a loga- 
rithmic-normal distribution. 

Figure  4 shows  the  distribution  at  every  floor  of 
the  maximum  pit  depth  detected  in  each  100  mm 
unit  length.  Though  the  maximum  pit  depths 
seemed  to  have  a slight  tendency  to  become  shal- 
lower at  upper  floors,  the  distribution  of  pit  depths 
was  regarded  as  being  uniform  throughout  the 
building.  The  actual  maximum  pit  depth  value  in 
the  present  survey  was  0.452  mm  which  was  de- 
tected at  6th  floor. 


Fig.  3.  Plot  on  a logarithmic-normal  distribution  diagram  of  all 
data  presented  in  Fig.  2. 

In  Figure  5,  the  maximum  pit  depth  data  for  unit 
lenghs  of  20  mm,  100  mm,  and  200  mm  were  evalu- 
ated from  extreme  value  analysis  on  the  basis  of 
Gumbel’s  double-exponential  distribution;  the  lin- 
earity of  each  plot  shows  that  this  distribution  ap- 
plies to  the  maximum  pit  depth  data  obtained  at 
the  unit  lengths  between  20  mm  and  200  mm. 

3.2  Minimum  Required  N and  .9 

In  the  practical  application  of  extreme  value 
statistics,  the  number  and  the  size  of  unit  samples 
for  the  pit  depths  survey  are  to  be  decided  prior  to 
the  destructive  (or  nondestructive)  inspection  of 
the  existing  structures. 

Reliability  of  the  extreme  value  prediction  de- 
pends on  the  following  three  equations; 

y =\nT  (3) 

r = y/5  (4) 

V{x ) = a\A  (N,n  )y^+B(N,n  )y+C (N,n )]  (5) 
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Maximum  pit  depth  (mm) 

Fig.  4.  Distribution  at  each  floor  in  the  building  of  the  maxi- 
mum pit  depth  detected  in  each  100  mm  unit  length. 


where  y = standardized  variable, 

T = return  period, 

5 = total  length, 
s =unit  length, 

V{x)  = variance  of  the  estimated  extreme 
value, 

A{N,n),  B{N,n),  C(A^,n)  = MVLUE 
coefficients, 
a = scale  parameter. 

It  is  obvious  from  these  equations  that  for  in- 
creased reliability,  hence  for  a minimized  V(x),  it  is 
required  to  increase  s (therefore  to  decrease  T) 
and/or  to  increase  N.  In  the  interest  of  economy, 
however,  there  must  be  a natural  limitation  to  the 
increase  of  N and  s for  the  extreme  value  survey. 

In  theory,  the  optimum  combinations  of  N and  .s 
can  be  determined  from  Eqs.  (3)  to  (5)  once  the 
distribution  parameter  (that  is,  a ratio  of  the  loca- 
tion parameter  to  the  scale  parameter,  a/A)  is 
reasonably  assumed  (that  is,  empirically  or  experi- 
mentally), and  the  extent  of  the  standard  deviation 
of  the  error  of  extreme  value  estimates,  a,  may  be 
expected  to  be: 

A =m(T  (6) 

where  A = location  parameter, 

o-=  standard  deviation  of  the  estimated 
extreme  value, 

m = assumed  number  (1,  2,  3,  etc.) 

That  is,  the  decision  may  be  made  to  reduce  a to 
1/m  of  A [6]. 


Maximum  pit  depth(mm) 
a.  Unit  length  20mm 


Maximum  pit  depth(mm) 
b.  Unit  length  100mm 


Maximum  pit  depth(mm) 
c.  Unit  length  200mm 


Fig.  5.  Double-exponential  distribution  type  of  extreme  value  probability  plots  for  i =20  mm,  100  mm,  and  200  mm. 
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Table  1 shows  the  result  of  extreme  value  analy- 
sis to  obtain  the  relationship  between  unit  length 
and  distribution  parameter.  In  the  analysis,  all  pit 
depth  data  obtained  from  the  whole  lengths  were 
brought  into  consideration.  It  can  be  seen  in  the 
table  that  the  distribution  parameter,  a/A,  de- 
creased with  increasing  unit  length,  ultimately  ap- 
proaching to  a definite  level  of  0.15. 


Table  1.  Location,  scale,  and  distribution  parameters  deter- 
mined from  all  maximum  pit  depth  data  obtained  in  the  survey 


Unit  length 
(mm) 

T 

N 

■^max 

(mm) 

A 

(mm) 

a 

(mm) 

al\ 

50 

177.6 

172 

0.575 

0.28415 

0.05621 

0.198 

too 

88.8 

86 

0.528 

0.30446 

0.04986 

0.16 

200 

44.4 

43 

0.506 

0.3221 

0.04866 

0.151 

Optimum  combinations  of  N and  T (hence,  s) 
required  for  controlling  a down  to  1/m  of  A were 
obtained  as  shown  in  Fig.  6 corresponding  to  a/A  of 
0.15  and  0.20.  It  stands  to  reason,  that  the  number 
of  unit  lengths,  N,  can  markedly  be  reduced  by 
loosening  the  requirement  for  the  reliability,  that 
is,  by  decreasing  m . N may  also  be  reduced  by  de- 
creasing a/A.  Though  N may  be  reduced  by  in- 
creasing s (hence  decreasing  T),  it  should  be  noted 
that  the  total  length  required  for  a pit  depths  sur- 
vey can  increase  on  the  contrary. 


Fig.  6.  Optimum  combinations  of  N and  T to  control  o-at  levels 
of  1/2,  1/3,  1/4,  or  1/5  of  A under  the  limitation  of  a/A  =0.2. 


3.3  Maximum  Pit  Depths  Prediction  Based  on  the 
Optimized  N and  s 

Figure  7 shows  the  results  of  the  maximum  pit 
depths  prediction  by  using  the  combinations  of  N 
and  s from  Fig.  6.  Unit  lengths  for  this  evaluation 
were  not  sampled  at  random  locations,  but  were 
ordered  from  one  end  of  the  lower  floors  upward. 
Thus  the  results  were  represented  by  correlating  to 
the  locations  wherefrom  those  unit  lengths  used  in 
the  prediction  were  removed.  It  can  be  seen  that 
the  scatter  of  estimated  extreme  values  decreased 
with  increasing  s.  As  would  be  expected,  the  scat- 
ter was  narrower  at  m = 5 as  compared  with  m -3. 

To  determine  the  general  tendency  in  the  reli- 
ability of  the  present  extreme  value  prediction. 
Table  2 was  developed  from  data  on  Fig.  7 by  tak- 
ing an  average  for  each  item.  The  results  indicate 
the  following: 

1)  The  assumptions  for  a/A  and  m,  that  were 
made  prior  to  the  analysis,  have  conservatively 
been  met. 

2)  Increased  m did  not  always  result  in  increased 
reliability  which  indicated  that  m-3  might  be 
reasonable  in  the  interest  of  economy. 

3)  Except  for  the  cases  of  5 = 10  mm,  the  maximum 
detected  pit  depth  value  of  0.452  mm  fell  closely 
between  the  average  of  estimates  ±lcr. 

4)  Although  the  level  of  s is  to  be  kept  as  low  as 
possible,  because  the  total  sampling  length  (ie., 
N xs)  may  increase  with  increasing  s,  at  least 
2.5%  of  the  whole  structure  is  to  be  subjected  to 
extreme  value  evaluation  at  a return  period  of 
500  or  the  less. 


4.  Concluding  Remarks 

Based  on  this  analysis,  the  following  conclusions 
can  be  drawn  concerning  the  optimum  conditions 
for  obtaining  a reasonable  extreme  value  predic- 
tion. 

1)  The  number  and  the  size  of  unit  samples  for  the 
extreme  values  survey  may  be  determined  so 
that  the  variance  of  the  extreme  value  estimates 
is  to  be  minimized  under  a definite  distribution 
parameter,  wherein  the  standard  deviation  of 
the  estimates  is  expected  to  be  1/3  of  the  mode 
of  distribution. 

2)  A return  period  of  500  with  a total  sampling 
length  amounting  to  2.5%  of  the  entire  parts 
may  be  desirable  for  an  increased  reliability  of 
extreme  value  prediction. 
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Fig.  7.  Distribution  at  each  floor  of  the  estimates  of  maximum  pit  depths  on  the  basis  of  the  optimum 
combinations  of  N and  T under  the  limitations  of  A =3cr  and  5cr. 


Table  2.  Summary  of  the  extreme  value  prediction 


Unit  length 

5 

(mm) 

Return 

period 

T 

Sample 

number 

N 

Location 

parameter 

Aav 

Scale 

parameter 

ttav 

Estimated 
maximum 
pit  depth 

^^max 

Aav 

^^max 

10 

888 

14 

0.241 

0.043 

0.534 

0.069 

3.50 

i ^ 

20 

444 

12 

0.269 

0.034 

0.478 

0.054 

4.98 

A “*  JU 

40 

222 

9 

0.288 

0.035 

0.467 

0.058 

4.98 

100 

88.8 

7 

0.313 

0.028 

0.428 

0.045 

6.94 

10 

888 

33 

0.218 

0.046 

0.559 

0.047 

4.63 

20 

444 

28 

0.266 

0.040 

0.491 

0.040 

6.67 

\=5(t 

40 

222 

23 

0.282 

0.039 

0.487 

0.039 

7.25 

100 

88.8 

18 

0.312 

0.034 

0.463 

0.032 

9.71 
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